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Abstract 

The analytical structure of some generalizations of the circle map 
is given. Also a generalization of off centre reflection is studied. The 
stability of Ito-Glass coupled map lattice is studied. 

Key words: Circle map, coupled map lattice, global Ito- Glass cou- 
pled map lattice. 

1 Introduction 

Circle map is one of the subjects which is both interesting mathematically 
and useful in applications [Glass et al 1989, Vandenmeer 2001]. It relates, in 
a beautiful way, some pure mathematical problems and applications. 

For example in [Glass et al 1989] it is shown that the simplest circle 
map f{6) = 9 + r mod 1, where r is a constant is related to a human 
heart disease (parasystol) where there are two types of heart pulses. This is 
represented by a symbolic sequence of 2-symbols (say 0,1). Such sequence 
defines another one called the reduced sequence where one starts with 1 and 
count the number of iterates till the next 1 appears and so on. This reduced 
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sequence has remarkable properties with some relevance to number theory: 
i) The reduced sequence contains, at most, three integers say p,m,n. ii) 
p = m + n. iii) m or n is odd. iv) Only one of p, m, n can succeed itself in 
the reduced sequence. These results agree with clinical observations. 

There are , however, some diseases e.g. tachyacardia which should be 
represented by coupled map lattice [Ito and Glass 1992]. Therefore in this 
paper we study some generalizations of the circle map analytically. The 
paper is organized as follows: In Section 2, the definition of the circle map 
is given, and some generalizations are introduced. Many generalizations Ito- 
Glass coupled map that is more suitable for modelling some heart diseases 
are introduced and studied in section 3. The conditions for the stability of 
its equilibrium solutions are derived. Some conclusions are summarized in 
section 4. 

2 Circle map and some generalizations 

In the following, some basic results about circle map will be reviewed and 
generalized. 

Definition 1: A circle map is a map / : S 1 — > S 1 , where S 1 is the unit 
circle. A map F : R — > R is the lift of the circle map / if 



where R is the real numbers and F o ix means the superposition of the two 
maps F and n, and 



7T O / = F O 7T, 



(1) 



7T : R -> S 1 , n(t) = exp(2rrit) V t G R. 



(2) 



In this case F (hence /) is of degree one if 




(3) 
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Definition 2: 

(a) The rotation number of a circie map / is defined by 

p(f, x) = lim 1 1 , F n = F o F n ~\ (4) 



(b) ff / is differentiable then the Lyapunov exponent of / is defined by 

1 Tl—l 

A = lim -£H/'(*;)I, (5) 

j=0 

where 

^+i = ffa), j = 0,1,2, .... (6) 



In general p(f, x), if it exists, depends on the initial point and forms a closed 
interval [Ito 1981]. 

The following is a known result yet for the sake of completeness we present 
a proof: 

Proposition 1: (a) If F(x) is monotonically increasing and continuous then 
p(f, x) is independent of x. 

(b) p(f q ) = qp(f) for all positive integers q. 

Proof. Since F is monotonically increasing and continuous (hence bounded) 
then lim ra ^ 00 (F n (a;) — x)/n exists. 

To show that it is independent of x assume, without loss of generality, 
< x < y < 1. Since F is monotonically increasing then there is a positive 
integer k such that 

F k (x) < F{y) < F k+ \x), 

hence 

F n+k -\x) < F n (y) < F n+k (x). 
Taking the limit as n — > oo, one gets 

lim F n (x) = lim F n (y). 
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Recalling that — 1 < x — y < 0, then 

Mm F "^- X = lim F ^-V. 

n— >oo n— >oo fl 

To prove (b) use 

F nq (x) - x , F nq (x) - x 
hm = q hm . 

n— >oo n gn— >oc qn 

This completes the proof. ■ 

Lemma 1 [Boyland 1986]: Let / : S 1 — > S* 1 be degree 1 continuous cir- 
cle map, then Vp, g positive integers, p/q E p(f, x) if and only if / has a p/q 
periodic point i.e. F q (x) = x + p. 

Definition 3: A circle map g is in class A if g is continuous, degree one 
and if its lift G is piecewise strictly monotone with precisely two critical 
points (g'(x) = 0) in (0, 1) one maximum at vri2 and one minimum at m\ and 
mi > m 2 . 

Definition 4: Let 

F = x+w+bp(x), p : R — > i?, p is continuous, = < w < 1, 

(7) 

such that 



(a) p(0)=p(l)=0. 

(b) If b > 1 then F is in class A. 

(c) If b < 1 then F is a homeomorphism. 

(d) is the restriction to the real line of a complex entire map. 

Theorem 1 [Boyland 1986]: Let F be as in Definition 4 and let / be the 
corresponding circle map then the boundaries of the regions in the b—w plane 
where p(f) is periodic (these regions are typically called Arnold tongues) are 
continuous curves. 
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Corollary 1: Let 



^2 



= x + w H sin(27ra;) + — sin(47r:r), (8) 

27r 47r 

then Arnold tongues for F(x) are continuous curves for r < 0.5. 
Proof. Comparing Eq. (7) and Eq. (8), then 

p(x) = sin(27ra;) + - sin(47rx), b = — . (9) 

2 27T 

It is direct to see that p(x) in Eq. (9) satisfies all the conditions for Definition 
4, hence using Theorem 1, the corollary is proved. ■ 

There are several ways to generalize the above results. One way is to use 
circle maps other than the standard sine circle map 

f(x) = x + w + 7^- sin(27rx), (10) 

27T 

as has been done in Eq. (8) and will be done in the following. Alternatively 
maps with degree other than one will be studied. This direction was pio- 
neered by Keener [Keener 1980] and extended by Glass and his group [Belair 
and Glass 1985, Bub and Glass 1995]. This direction will be discussed further 
later on. 

Definition 5 [Au 2001]: The off centre reflection map / : S 1 — > S* 1 is 
defined by 

/ = x + w - 2{tan~ 1 [sin ] - x}, < r < 1, < w < 1. (11) 

cos 2ttx — r 



The following proposition was briefly proposed in [Au 2001] paper. Here 
the full proof is given: 

Proposition 2: The off centre reflection map f(x) in Eq. (11) can be 
written as 

OO 

f(x) = x + w -2^~2 1 —sm(2irkx). (12) 
fc=i 



Proof. The third term of f(x) in Eq. (11) is an odd function hence it 
is representable as Fourier sine series hence 

oo 

f(x) = x + w — 2 b k (r) sin(27r£;x), 
k=i 

where 

db k 2 r cos(k — l)u — cos(k + l)-u 

-aw. 



(9r 7T y 1 — 2r cos -u + r 2 
Define 

2 r cos(fcw) 

*k — — z ^ : — n c 

if J 1 — 2r cosm + r z 

and use 

1 



(1 — rexp(iu)) x (l — r exp(— iu)) 1 



1 — 2r cosu + r 2 
Expanding in powers of r, one gets 



OO 

„2n+k _ T 



n=0 



thus one gets 



which completes the proof. 



Corollary 2: The function in Eq. (12) is phase of the restriction on the 
unit circle of the complex function 

z 2 (l - rz) 

h{z) = exp{tw)^ LfZ, 0<r<l. (13) 

z — r 



The proof is by taking the logarithm of both sides and expanding in 
powers of r. 

It is straightforward to derive the following result [Au 2001]: 
Proposition 3: 
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(a) For < r < 1/3, the off centre reflection map / in Eq. (12) is mono- 
tonic hence its rotation number is a single point independent of the 
initial value. 

(b) There is a saddle node bifurcation at 

= — , cosa r = r. (14) 

7T 

(c) The off centre reflection map /, Eq. (12) has a 2-cycle if w = or 
w — 1/2 and a period doubling bifurcation at r = l/y/h to a 4-cycle. 

Now we consider maps of degree M(M ^ 1) i.e. 

f(x + l)=f(x)+M. (15) 




Definition 6 [Belair and Glass 1985]: The rotation number for a degree 
M map is defined by 

- rt— 1 

p{f,x)= lim _^[F(^)-^], (16) 

j=0 

where F is the lift of /. 

Proposition 4: [Belair and Glass 1989] 

(a) Let 

f(x, w) = g(x, w) + w, (17) 

be a circle map with degree M then Vw real number there is at least 
(Sj=o M^ -1 ) different values between w, w + 1 for which x is a periodic 
point of f n . 

(b) If g is strictly monotonic, g = g(x,b), b is a real parameter, then there 
is a continuous curve w(b) along which a; is a periodic point with period 
n. 

Definition 7 [Keener 1980]: Let F be the lift of a circle map with degree 
M. Following Keener [Keener 1980] the following properties are defined: 
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F satisfies property FI if F(x) is piecewise continuously different iable 
function defined on the interval [0,1] with a single discontinuity at x — 

9 e(o,i). 

F satisfies property FII if it satisfies FI and 

dF 

x^6 =>• < — < oo. (18) 
ax 

F satisfies property Fill if it satisfies FII and 

lim F(x) = 1, lim F(x) = , F{6) = or 1 9 e (0, 1). (19) 

where F has a single discontinuity at x = 9 e (0, 1). 

F satisfies property FIV if it satisfies Fill and F(0) > F(l) (non- 
overlapping case). 

F satisfies property FV if it satisfies Fill and F(l) > F(0) (overlapping 
case) . 

Definition 8: The system (6) is chaotic in the standard sense if its Lya- 
punov exponent, Eq. (5) is positive. It is chaotic in the sense of Keener if its 
rotation number, Eq. (4) depends on x. 

In the following we present two results. The first is about the rotation 
number of a generalization of the truncated off centre reflection, Eq. (12). 
The second is a proof that the two definitions of chaos in Definition 8 are 
inequivalent. Define 

F(x) = Mx + w- (^-) sin 2irx - (-^-) sin(47r:r). (20) 



Proposition 5: If the following conditions are satisfied then the rotation 
number of F(x) in Eq. (20) is 1/2 

M + W > , + (£) + (£),„ + (£) + (£)<!, 

a + b < M<1, (21) 
w > 9, M + w <l + 9. 



8 



Proof. Equation (20) implies 



Mx + w H h — > F(x) > Mx + w 

2-k A-k 



a b 
2tt ~ 4tt' 



thus 



1 + (a/2vr) + (6/4?r) - w 
M 



> > 



1 - (a/27r) 



M 



(b/4ir) - w 



where 9 is denned in Eq. (19). The first two conditions in Eq. (21) guarantee 
that < 9 < 1. 

The condition a + b < M implies that F(x) satisfies FII. The condition 
M < 1 implies that F(0) > F(l). This guarantees that the discontinuity of 
F is unique. Thus F satisfies FI-FIV hence by Keener's Lemma 3.1 [Keener 
1980] the rotation number of / is a single point. 

To prove that the rotation number equal to 1/2 use the third set of 
conditions in Eq. (21) which implies 



hence [Holmgren 1996] F 2 has a fixed point in (0, 9) i.e. the rotation number 
equal to 1/2. ■ 

To see that the conditions in Eq. (21) are consistent let 



Proposition 6: The two definitions of chaos (c.f. Definition 8) in the sense 
of Keener and the standard sense are inequivalent. 

Proof. Consider the map (20) with the following conditions 



F(0, 0) C (9, 1), F(9, 1) C (0, 9) =► F 2 (0, 9) c (0, 9), 



M = 0.8, w = 0.6, 9 = 0.5, a = b = 0.1. 



Now the following proposition can be proved. Let 



F = Mx + w mod 1. 



(22) 



2 - w > M > 1, < w < 1. 



(23) 
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The first condition in Eq. (23) implies that F satisfies FI, FII, Fill, FV 
(overlapping condition F(l) > F(0)). The second condition implies that x = 
(1 — iu)/(M — 1) is a fixed point and the third one implies < 9 = (l — w)/M. 
Since x — (1 — w)/(M — 1) is a fixed point of F it belongs to 1. Using Eq. 
(23), one gets 

F(l) = M + w-l<^-^, 

then I is not subset of [F(0), F(l)]. Moreover since x = is not a fixed point 
of F, then I ^ [0, 1]. Thus by Theorem 1 the rotation number of / is a single 
point. Yet Eq. (22) implies that Lyapunov exponent of F is positive. This 
completes the proof. ■ 



3 Ito-Glass coupled map lattice 

Ito and Glass [Ito and Glass 1992] have proposed a coupled map lattice which 
is more suitable for modelling the heart. The equations are given by 

i— 1 n 

< +1 = -9(4) + *[£ + E f( x % ( 24 ) 

j=i j=i 

where i = 1, 2, n. The homogeneous equilibrium solution is given by 

x = -g(x) +nlf(x). (25) 

Define a and c by 

a = l^c=^l. (26) 
dx ax 

Using inhomogeneous analysis they were able to show that the characteristic 
polynomial for the equilibrium solution Eq. (25) is 

(-1)»{A(A + c)« - [(1 + a)X-a + c] n } = Q 



A - 1 

We begin by studying the simplified system: 

i— 1 n 



f = ( 28 ) 



k=l k=i 
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The homogeneous equilibrium solution is given by 

= nf(6). (29) 
Using Eq. (27) the characteristic polynomial for the solution (29) is 



A - 1 



0, 



where a = f'(0), 9 is given by Eq. (29). Using l'Hopital rule it can be 
written as 

A™ - [n/(n + 1)](1 + a)[(l + a)A - af' 1 = (30) 

Since we are interested in modelling the heart we are interested in the 
case n is large. Thus we have: 

Proposition 7: 

1. The equilibrium solution, Eq. (29) is unstable if any one of the following 
conditions are satisfied: 

(i) a > 0, (1 + a) n+1 > a(n + l)/[n(l + a)]. 

(ii) a < —2, n is sufficiently large. (31) 

(iii) — 2 < a < — 1, nis odd and sufficiently large. 

2. There are no real roots for Eq. (30) satisfying |A| > 1 if — 1 < a < 0, 
and n is sufficiently large. 



3. Define 



(-1)*+Xn-1)! , ^ n _ k 



a = 1, a k+1 = {n+1)kKn _ k _ lV a (1 + a)- , fc = 0,1,2, ...,n-l, 

(32) 

then if —1 < a < and n is sufficiently large such that a k +i < 1 Vfc = 
0, 1, 2, n — 1, then the equilibrium solution (29) is stable. 

Proof. We will prove (i) by showing the existence of a root A > 1 for Eq. 
(30). If there is a root A»a/(a+ 1), then Eq. (30) can be approximated by 

77 77 

A n ~ (1 + a) n \ n - 1 A ~ (1 + a) n > 1 for n > 1. 

n+V n + V 
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The same argument is valid for a < — 2 and n sufficiently large hence (ii) is 
proved. 

To prove (iii) we will show that there is a positive integer k > 1 such that 
the LHS of Eq. (30) is positive at A = —k + 1 and negative at A = —k. Since 
n is assumed odd the first term in Eq. (30) is always negative for A < while 
the second term is positive (notice that 1 + a < 0). 

For n sufficiently large the positivity at A = — k + 1 implies 

k \a\ > 2(k - 1), 

while the negativity at A = —k implies [k + 1) \a\ < Ik. Thus whenever 
2(k -l)/k< \a\ < 2k/ (k + 1) there is a root A e [-k, -k + 1]. 

To prove the second part of the proposition assume — 1 < a < 0. The 
proof is by direct substitution for a = 0, a = — 1. For these two values the 
homogeneous solution (29) is stable. 

For — 1 < a < it is direct to see that, in this case, Eq. (30) does not 
admit a negative solution. Thus the only remaining possibility is a A > 1 
solution. To show that this is not possible, define /(A) to be the right hand 
side of Eq. (30) and f k (\) to be the k-th derivative of / with respect to A. 
Then it is direct to see that 



f k (l) = n(n-l)(n-2)...(n-k+l) 



{n-k){l+af 
n+ 1 



fc+i 



> VA; = 0,1,2, 



Consequently, using Taylor series we have /(A) > VA > 1. Therefore there 
are no real roots A > 1 for Eq. (30) in the case — 1 < a < 0. 

We prove the third part of the proposition using Jury test [Murray 2002] 
so we define 

b n „ k = a k — a n a n - k , < k < n - 1, 
c n -k = KK-k - hh+i, < k < n - 2, 
d n - k = c n c n -. k - c 2 c k+2 , < k < n - 3, 

g n __ fe = VnVn-k - p n -3Pk+n-3, < k < 2. 

The conditions we assume in part (3) of the proposition implies 



bn—k — C^kj Cn—k — &ki d n — k — Oi k) Qn—k — (%k> 
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hence 

\a n \ < 1, < \b n \ , |c 2 | < |c n | , |d 3 | < , |g„_ 2 | < |g n | • (33) 

Since 

/(D - i-^U/(i)>o, 

(-l) n /(-l) = 1 + ( 2a + I)"" 1 > 0, V - 1< a < 0, 

then the conditions of Jury test are satisfied. ■ 
Applying Proposition 7 to the circle map 

f(x) = x + w + ■^-ean(2irx), (34) 

2,71 

one gets that if k < 1 then the homogeneous solution for the Ito-Glass CML 
for the circle map (34) is unstable if (31) is valid. 

We introduce the following local version of Ito-Glass coupled map lattice 
(IGCML) 

xf 1 = (1-D)f(x\) +j [f{xt\) + /(4k)] , i = 1, 2, n and < D < 1. 

(35) 

The homogeneous equilibrium solution is given by 

x = /(*). (36) 

Proposition 8: The equilibrium solution (36) is stable if \f'(x)\ < I. 
Proof. Assume 

x\ =x + el |e?| < r]\/i = l,2,...,n. 
Using the conditions < D < 1, \f'(x)\ < 1 then 

e\ = (1 - D)f'(x)e° 1 + ^f'(x)e° 2 =>\e\\<r,, 
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hence 



e 1 2 = (l-D)f'(x)4 + ^f(x)(4 + e\) =► |4| < rj : 



and so on. Thus < rj \/i = 1,2, ...,n and Vt > 0. This completes the 
proof. ■ 

We introduce another form of global Ito-Glass coupled map lattice which 
is significantly closer to those of Kuramoto, Winfree etc.. [Ahmed and Hegazi 
2002] is 



x. 



t+i 



k 



n 



i-l 



E/(4 +1 ) + £/(*$) 

j=l j=i 

< wt < 1 Vi= l,2,...,n, 1 > A; > 0. 



(37) 



The equilibrium solution is given by 

w + kf(x) = 0. 



(38) 



Definition 9: A dynamical system is persistent if x° > Vi = 1,2, ...,n 
implies lim^oo inf (x-) > ,Vi = l,2,...,n. 

Proposition 9: 

(a) For sufficiently large n, the solution (38) is stable if f'(x) < 0. 

(b) The system (37) is persistent if f(0) > 0. 

Proof. Applying l'Hopital rule to Eq. (27), one gets that the character- 
istic polynomial for the homogeneous solution of the general Ito-Glass CML 
is 

(A + c) n + n\(X + cf- 1 - n{\ + a)[(l + a)X - a + cf- 1 = 0. (39) 
The system (37) corresponds to 

k 

c = — 1, a = — fix). 
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Substituting in Eq. (39), one gets A = 1 or 
A 



1 n 
+ 



n+1 n + 1 



1 + -f(x) 
n 



for n sufficiently large 

A ~ exp(fc/'(a;)), 

which proves (a). 

To prove (b): Assume 1 x\ > for some t and Vi = 1, 2, n and use 
Eq. (37) hence 



x 



t+i 



xl + w i + f(0)(-) 
n 



i—l n 

.j=i 



3=1 



x\ +1 ~ x*+wi + / / (0)(-)4 >xi, 



4+ 1 ~ 4 + u , 2 + /'(o)( 



n 



x\ +1 + 



> x 2 , 



and so on. This completes the proof of the proposition. 
Another form of Ito-Glass coupled map lattice is 



*r = (i-k)f(xD + 



k 



n — 1 



i-l 



(40) 



< Wi < 1 Vi = l,2,...,n , 1 > A; > 0. 



The homogeneous equilibrium solution is given by rr = /(x). If 1 then 
it is direct to see that the equilibrium solution is stable if |/'(a:)| < 1. For 
general k (1 > k > 0) and sufficiently large n one has 

Proposition 10: For general k (1 > k > 0) and sufficiently large n the 
roots of the characteristic polynomial of Eq. (40) are 



A = — c, and A ~ exp 



-kf'(x) 



1-A 

A + c 



, c=(-l + fc)/ / (x). (41) 



Moreover, if —1 < f'(x) < then the equilibrium solution x = f(x) is stable. 
Proof. The proof of Eq. (41) is similar to that of part (a) in Proposition 9. 
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Now assume that —1 < f'{x) < 0, then |c| < 1. Thus instability may 
arise only from the transcendental equation in Eq. (41). If A > 1, then 
—kf'(x)(l — A) /(A + c) is negative hence exp[— kf'(x)(l — A)/(A + c)] < 1, 
contradiction. Since exp[— kf'(x)(l — A)/(A + c)] > 0, then A < — 1 is not 
possible. ■ 

4 Conclusions 

Some generalizations of the circle map are given. The two definitions of chaos 
(in the sense of Keener and the standard one) are found to be inequivalent. 
Ito-Glass coupled map lattice is more suitable for modelling the heart. Many 
forms of Ito-Glass coupled map lattice are introduced and studied. The 
conditions for the stability of its equilibrium solutions are derived. 
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